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Abstract — This paper presents a novel framework for modeling 
the uplink intercell interference (ICI) in a multiuser cellular 
network. The proposed framework assists in quantifying the 
impact of various fading channel models and state-of-the-art 
scheduling schemes on the uplink ICI. Firstly, we derive a semi- 
analytical expression for the distribution of the location of the 
scheduled user in a given cell considering a wide range of 
scheduling schemes. Based on this, we derive the distribution and 
moment generating function (MGF) of the uplink ICI considering 
a single interfering cell. Consequently, we determine the MGF 
of the cumulative ICI observed from all interfering cells and 
derive explicit MGF expressions for three typical fading models. 
Finally, we utilize the obtained expressions to evaluate important 
network performance metrics such as the outage probability, 
ergodic capacity, and average fairness numerically. Monte-Carlo 
simulation results are provided to demonstrate the efficacy of the 
derived analytical expressions. 

I. Introduction 

Explosive growth in the demand of high quality wireless 
data services compel the network designers to utilize spectrum 
more aggressively which on one side enhances the spectrum 
efficiency, whereas on the other side it enhances the intercell 
interference (ICI) which is an alarming bottleneck in the 
telecommunication growth paradigm. The allocation of the 
same frequency bands across neighboring cells produces inde- 
terministic ICI which is highly dependent on the statistics of 
the channel characteristics and on the dynamics of the schedul- 
ing decisions. In this context, it is of immense importance for 
the system designers to accurately characterize and investigate 
the behavior of the ICI which helps in gaining more theoretical 
insights, quantifying various network performance metrics 
and developing efficient resource allocation and interference 
mitigation schemes. 

Orthogonal frequency division multiple access (OFDMA) 
has been recently adopted as the multiple access scheme for 
the state-of-the-art LTE and WiMAX cellular technologies. 
In OFDMA, a wide-band frequency-selective fading channel 
is decomposed into a set of orthogonal narrow-band sub- 
carriers. The orthogonality among the subcarriers per cell 
makes the intra-cell interference almost negligible. However, 
with universal frequency reuse among cells (i.e., all cells 
use the same set of subcarriers), the ICI at each subcarrier 
may cause severe degradation in the network performance. 
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In OFDMA networks, the subcarriers are allocated adaptively 
among users per cell based on a predefined scheduling scheme. 
Moreover, each subcarrier is allocated to only one user per cell 
assuming BSs are equipped with single antenna and, thus, the 
number of interfering users on each subcarrier is rather limited. 
Therefore, the cumulative ICI on a given subcarrier may not 
be modeled accurately as a Gaussian random variable (RV) by 
simply invoking the central limit theorem. 

Several recent studies considered the modeling of ICI in 
the downlink where the location of interferers is typically 
deterministic. A semi-analytical distribution for the signal-to- 
interference-noise ratio (SINR) has been derived in Q] under 
path loss and log normal shadowing for randomly located 
femtocell networks. In 0, the applicability of the Gaussian 
and binomial distributions for modeling the downlink ICI 
is investigated. In 0, the optimal threshold is derived for 
fractional frequency reuse (FFR) systems assuming ICI as 
Gaussian RV. In 0, the authors derived the distribution of 
the ICI under log-normal shadowing and Rayleigh fading. 
The distribution of ICI is shown to highly deviate from the 
Gaussian distribution in OFDMA networks. 

In comparison to downlink, the nature of uplink ICI is 
different in various aspects that include the following: (i) Due 
to the implicit symmetry and fixed locations of the BSs in 
the typical grid-based downlink network models, the number 
of significantly contributing interferers typically remains the 
same irrespective of the position of the mobile receiver. Also, 
it has been shown in [2\ that the strongest interference is 
generated by two closest interfering BSs irrespective of the 
mobile receiver location. However, the number of significantly 
contributing interferers in the uplink cannot be quantified 
at a given instant due to the highly varying locations of 
the interfering mobile transmitters; (ii) Conditioned on the 
location of the desired mobile receiver within a cell, the 
exact distance of the interfering BSs can be calculated in 
the typical grid-based downlink network models. However, 
knowing the location of the BS receiver in the uplink does 
not help in determining the exact location of the interfering 
mobile users; (iii) In the uplink, cell edge and cell center 
mobile users are subject to the same amount of interference 
on a given subcarrier, which is the interference received at the 
BS. Whereas the same is not true for the downlink in which 
cell edge users experience higher interference coming from 
the nearby BSs 0, 0. 

Some worth mentioning research works for the uplink ap- 
pear in 0, 0-1(9]. In 0, the authors developed an analytical 
model for subcarrier collisions as a function of the cell load 
and frequency reuse pattern. They derived an expression for 
the SINR in the uplink and downlink, ignoring the effect 
of shadowing and fading. In 0, the authors developed an 
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analytical expression for the subcarrier collision probability 
considering non-coordinated schedulers. In (8), the authors 
modeled uplink ICI in an OFDMA network as a function of 
the reuse partitioning radius and traffic load assuming arbitrary 
scheduling. In (9). the authors presented a semi-analytical 
method to approximate the distribution of the uplink ICI 
through numerical simulations without considering the impact 
of scheduling schemes. 

In this paper, we propose a novel theoretical framework to 
derive the statistics of the uplink ICI on a given subcarrier 
as a function of both the channel statistics (i.e., path loss, 
shadowing and fading) and scheduling decisions. The pro- 
posed framework can be also extended to typical downlink 
scenarios as explained in flTOl . The framework is generic in the 
sense that the derivations hold for generalized fading channels 
and various scheduling algorithms. We start by deriving the 
distribution of the location of the scheduled user in a given 
cell. We then derive the distribution and moment generating 
function (MGF) of the ICI considering a single interfering cell. 
Next, we derive the MGF expression for the cumulative ICI 
experienced from all interfering cells over generalized fading 
channels, and present explicit expressions for three practical 
fading models. Finally, we demonstrate the importance of the 
derived expressions by utilizing them to evaluate important 
network performance metrics such as outage probability and 
ergodic capacity. 

The remainder of this paper is organized as follows. Sec- 
tion II presents the system model and the main steps of the 
proposed framework. In Section III, the distribution of the 
scheduled user location is derived for different scheduling al- 
gorithms. In Section IV, the distribution of the uplink ICI from 
one neighboring cell is derived. The MGF of the cumulative 
ICI from all interfering cells is determined in Section V and 
utilized in Section VI to evaluate three network performance 
metrics. Finally, numerical and simulation results are presented 
and analyzed in Section VII, and conclusions are drawn in 
Section VIII. 

Notation: Exp(A) represents an exponential distribution with 
parameter A, Gamma(m s , m c ) represents a Gamma distri- 
bution with shape parameter m s and scale parameter m c . 
t^G( m c, ms, represents the Generalized-/C distribution with 
fading parameter m c , shadowing parameter m s and average 
power 0. r(.) represents the Gamma function. P(A) denotes 
the probability of event A. /(.) and F(.) denotes the proba- 
bility distribution function (PDF) and cumulative distribution 
function (CDF), respectively, [a, b] denotes a discrete set of 
elements which ranges from a to b. Finally, E[.] denotes the 
expectation operator. 

II. System Model and Proposed Framework 
A. Description of the System Model 

We consider a given cell surrounded by L interfering 
neighboring cells. For analytical convenience, the cells are 
assumed to be circular with radius R. Each cell I is assumed 
to have U uniformly distributed users. The frequency reuse 
factor is assumed to be unity with each subcarrier reused in 
all cells. The bandwidth of a subcarrier is assumed to be less 
than the channel coherence bandwidth, thus, each subcarrier 
experiences flat fading. Time is divided into time slots of 
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Fig. 1 . Geometrical illustration of dividing the cellular network into multiple 
rings of non-uniform width Afe. 

length smaller than the channel coherence time and, thus, the 
channel variation within a given time slot is negligible. 

Generally, the scheduling strategies can be broadly catego- 
rized into two classes; (i) rate maximization (i.e., rate adap- 
tation) while transmitting with constant/maximum power; (ii) 
power minimization (i.e., power adaptation) while achieving a 
fixed data rate. In this work, we focus on rate adaptive schemes 
where users transmit with their maximum power in order to 
maximize their rate depending on the existing channel and 
interference conditions. Therefore, for the scope of this paper, 
we assume that all users transmit with their maximum power 
F mm on a given subcarrier with rate adaptation depending on 
their channel qualities. At this point, it is however important to 
emphasize that this is not a limitation and the approach can be 
extended for uplink power control mechanisms as discussed 
in iflOl . The instantaneous signal to noise ratio (SNR) 7 of 
each user can then be written as follows: 



7 = PmexC- 



(1) 



'-, C is the path loss constant, r is the user 



where K = - 

distance from its serving BS, ip and r} denotes the shadowing 
and small scale fading coefficient between user and BS on a 
given subcarrier, respectively, j3 is the path loss exponent, a 
denotes the thermal noise at the receiver and £ is the composite 
fading. Note that all users are assumed to be associated with 
their closest BS 0, CD, therefore r < R. 

Each cell is divided into K concentric circular regions. 
Since path loss decays exponentially from cell center to cell 
edge, therefore, we consider discretization of cellular region 
in such a way that the path loss decay within each circular 
region remains constant or uniform. The main motivation for 
dividing the cell into a discrete set of circular regions relies 
on the fact that the channel statistics of the users located in a 
given circular region become relatively similar especially for 
large values of K. More explicitly, the characterization of the 
circular regions can be demonstrated as follows: 



logio^fc = 



k + 10/31ogirjrfc-i 
10/3 



r <R, 



(2) 



where k is the path loss decay within each circular region 
[dB]. Due to the exponential nature of the path loss, it varies 
rapidly near the cell center than at the cell edge, therefore, 
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(i) each of the k th circular region bounded by two adjacent 
rings, i.e., r k and r k -\ possess non-uniform width = ?*& — 
Tk-i', (ii) the number of circular regions are high in the cell 
center than at the cell edge; (iii) the average number of users 
located within k th circular region bounded by ring r k and 
rfc_i are considered to be located at Note that, this is an 
approximation which is required for deriving the analytically 
tractable model of ICI and in any case it is not required for 
the Monte-Carlo simulations. The average number of users in 
each ring k (for analysis) can then be given as: 



U(r 



' k-i) 



R 2 



k = l,2,---,K, 



(3) 



It is important to note that Uk can be a fraction; therefore, we 
round off the fractional part of users. 

B. Main Steps of the Proposed Framework 

In order to characterize the statistics of the uplink ICI for 
generalized fading channels and various scheduling schemes, 
the proposed framework mandates the following steps: 

i) Derive the distribution / rsel (r) of the distance of the 
allocated user r sc i in a given cell from its serving BS 
based on the deployed scheduling scheme. 

ii) Derive the distribution ff aBl {f) of the distance between 
the allocated user in a neighboring interfering cell and 
the BS of the cell of interest f sc \. 

iii) Derive the distribution fx t (x) of the interference from 
the allocated user in neighboring cell I to the BS of 
interest, i.e., Xi = Kr~j^x, where \ denotes the com- 
bined shadowing and fading, i.e., composite fading of the 
interference. Finally, derive the MGF of the cumulative 
ICI, i.e., Y = YjI=i Xi from the scheduled users on a 
given subcarrier in all L interfering cells. 

III. Distribution of the Scheduled User Location 

Considering the dependence of the uplink ICI on the 
location of the scheduled users in the neighbor interfering 
cells which in turn depends on the deployed scheduling 
schemes, we derive in this section the distribution of the 
distance between the scheduled user and its serving BS in 
a given cell (i.e., the probability mass function (PMF) of 
r so i) considering the following five scheduling algorithms: 
greedy scheduling, proportional fair scheduling, round robin 
scheduling, location based round robin scheduling, and greedy 
round robin scheduling. 



A. Greedy Scheduling Scheme 

Greedy scheduling is an opportunistic scheme that aims at 
maximizing the network throughput by taking full advantage 
of multiuser diversity. However, it suffers from low fairness 
among users which makes it less attractive for network 
operators. The procedure for determining the PMF of r so \ 
considering greedy scheduling is divided into two steps: 
Step 1 (Selecting the user with the highest SNR in ring 
k): Since the path loss decay within each circular region 
is considered to be uniform, we approximate the distance 
of all users located within fc th circular region by ring r k 
for analytical tractability as we already mentioned in Sec- 
tion II. In this step, we select a user with maximum SNR 



in each ring k which possess Uk users. Thus, selecting a 
user in a ring k is equivalent to selecting the user with 
maximum channel gain among all the users in ring k, i.e., 
( k = max{Ci, <2, • • • Xi,'-' > Cu k }, where Q is the composite 
fading channel gain between user i and its BS on a given 
subcarrier. The CDF and PDF of the maximum channel gain 
C,k can be written as follows: 



(4) 



»=i 



fc (a) = E fa n *c< (c^ = «*/c(co mtk)) 

(5) 

To consider path loss, we now perform a transformation of 
RVs using (Q3, 7fc = KrZ^ Xk> where, j/. is the selected user 
SNR in each ring k. The CDF and PDF of can then be 
written as follows: 



u k -l 



^ fc (7 fc ) = n F C.( /? " 1 ^) i ' 4d ( F c(7 fe rfA- 1 )) , (6) 

i=l 

-, "ft "fc 



Kr h 

k j=l 
i.i.d u k 



Kr, 



Uk—l 



(7) 



Step 2 {Selecting the user with maximum SNR among K 
rings): In this step, we compute the probability of selecting 
the A: th ring among all other rings. It is important to note that 
this is equivalent to selecting the ring k which possesses the 
user with the highest SNR among all rings. Conditioning on 
7fc, the PDF of r sc i can be written explicitly as follows: 

K K 

P(r sci = r k \ lk ) = p{ ll < lk )= 1] f 7,(7fc), 

(8) 

By averaging over the distribution of 7/-, the final expression 
for the PMF of r sc i is 

P(r scl = r k )= [ ( f[ F 7i (7fe) J f lk (^k)dj k , (9) 

Using (O, (O can be written for i.i.d. case as follows: 

P(r se i = r k )= r fl (f^K- 1 ))^ x 

J ■ 1 - / 1 



(10) 



Kr-J 



(F c ( 7fc rf^- 1 )) k 



where r sc \ £ [0, R). The results in (fTOb are generalized for 
any shadowing and fading statistics. Even though ( [Tol l is not 
a closed form expression, the integration can be solved ac- 
curately and efficiently using standard mathematical software 
packages such as MAPLE and MATHEMAT I CA. 
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B. Proportional Fair Scheduling Scheme 

The proportional fair scheduling scheme allocates the sub- 
carrier to the user with the largest normalized SNR (7/7) lfl2l . 
where 7 and 7 denote the instantaneous SNR and the short 
term achieved average SNR of a given user, respectively. In 
other words, the selection criterion is based on selecting a 
user who has maximum instantaneous SNR relative to its own 
average SNR. The distribution of r sc \ can be derived as: 
Step 1 (Selecting the user with maximum normalized SNR 
in ring k): In this step, the performance of proportional fair 
scheduling scheme is independent of the path loss factor if 
users are moving relatively slowly, i.e., their path loss remains 
nearly the same on a short term basis. In this case, the problem 
of selecting the maximum normalized SNR in a ring k can be 
written as: 



1 Ci C2 

Cfe = max { j-, 

. Ci C2 



^ , (ID 



Ci Cuk 
Ci Cu k 

where Q = J °° Cif& (CO^Ci I s the average of the com- 
posite fading channel and Cfe is the maximum normalized 
composite fading channel gain in ring k. For i.i.d. aver- 
age composite fading gains of the users located in ring k, 
i.e., C = Ci = C2 = ■■■C,u k , the problem of selecting 
the user with maximum normalized channel gain reduces 
to selecting the user with the maximum channel gain, i.e., 
Cfe = max{Ci, C2, ■ • ■ , Ci, ■ ■ ■ Xu k }- Thus, for any ring fc, the 
CDF and PDF of the selected SNR 7 fe = Kr^Ck can be 
written as: 



(12) 



; = 1 



/3 u k 



i=^,i^3 



(13) 

The short term average SNR of the selected user in ring fc, i.e., 
7 fc can then be computed as % = J °° 7fe/ 7fe (7fe)c?7fe. Finally, 
the normalized selected SNR in each ring k can be defined as 
Cfe = ^ an d performing a transformation of RVs, the CDF 
and PDF of Cfe can be given as i*£ k (£fc) = Fjkilk^k); and 

= IkfikilkCk), respectively. 
Step 2 {Selecting the ring k with maximum normalized SNR 
from the K rings): Once we characterize the PDF and CDF 
of Cfe, the probability of selecting any ring k can be written 
using ([8]) as follows: 



P(r se \ = r k ) 



Hi n fuaol/fo 



(CfeRfe. (14) 



C. Round Robin Scheduling Scheme 

Round robin scheduling is a non-opportunistic scheme 
where a user is selected randomly within a time slot. As 
each user has equal probability of allocation, round robin is 
a strictly fair scheduling scheme. The round robin scheme 
provides maximum fairness among users and may serve as a 
lower bound in terms of network throughput which is useful 



in calibrating the performance of other scheduling schemes, 
however, the resulting network throughput is significantly low 
which makes it less attractive for practical implementations. 
The PMF of the scheduled user location can then be given as 



P(r se i 



r k = 



U 



(15) 



D. Location Based Round Robin Scheduling Scheme 

Location based round robin is another non-opportunistic 
scheduling scheme which do not require any channel state 
information, however, it requires the location information 
of the users. Even though location based scheduling is not 
common in practice, the location of each mobile user can be 
determined at the BS using global positioning system (GPS) 
or estimate based on a power measurement of pilot signals 
from the surrounding beacons, e.g., using triangulation based 
techniques. In this regard, there are variety of techniques 
available in the literature which demonstrate how the location 
of users can be evaluated at the BS (see JT3), lfT4l and the 
references therein). Moreover, the users located in different 
circular regions can also be classified based on the long term 
average SNRs, i.e., by computing SNR thresholds for different 
distances (rings) O which is a common technique in FFR 
systems to distinguish between cell-edge and cell center users. 

In this scheme, we consider W time slots during which 
the distance of the users from their serving BS will remain 
approximately the same. For simplicity, the number of time 
slots W is set equal to K. At a given time slot T w , we select 
any arbitrary user from a specific ring (analysis) and circular 
region (simulations) starting from the cell center. We continue 
to allocate the users by accessing the circular rings sequentially 
from cell center to cell edge. At this point, it is important to 
emphasize that all cells are considered to be time synchronous 
in allocating the users from particular rings, i.e., at a given time 
slot T w all cells are selecting the w th ring. Thus, the PMF of 
r so i for a given time slot w denoted by P(r sc \ = r£ w ) can be 
given as: 



P(r, 



scl 



') 




(16) 



E. Greedy Round Robin Scheduling Scheme 

Greedy round robin is an opportunistic scheduling scheme 
which captures the multiuser diversity while maintaining some 
fairness among users. In this scheme, we consider W = K 
time slots during which the distance of the users from their 
serving BSs remain nearly the same, however, the small scale 
fading gain on the considered subcarrier may vary from one 
time slot to the other. We select the user with maximum SNR 
in each time slot T w , however once a user is selected from a 
ring, all users located in that ring will not be scheduled for 
transmission for the next K — 1 time slots. Note that all BSs 
are considered to be time synchronized in terms of scheduling. 
Clearly, the probability of allocating a ring k at T% can simply 
be given by (TTOb . However, the probability of selecting a ring 
k at T2 is a dependent event and can be derived using Bayes 
theorem as follows: 



K 



P(r s , 



J2 P 



r se i = r k 



r sc \ 



.Ti 



P(r sel =rf), (17) 
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where, 



P 



( r se \ = r k T2 \ 




\r se l = Tj T l ) 





K 



i^k 



(18) 



UkfdlkrlK 2 ) 



Fdlkr'lK- 1 ) 



d~f k 



Since the probability of allocating any ring k within time slot 
T w depends on all previous states, therefore, the principle of 
Markov chain transition probabilities is not directly applicable. 
For more clarity, the probability of selecting a ring k at T3 is 
given as follows: 



K K 



r sc \ = r k 



r 2 



p 



r sc \ 



r sc i 



P(r s 



r Tl ) 
(19) 



where, P 1 



- T 3 



- Tl J is given at the top of next 



r S cl=rj T2 nr B , 

page in (120) . In general, the probability of selecting any ring 



fc at a time slot T w , i.e., P(r s , 



can be written as: 



K K 



EE" " E ^el = r-)P 
_/ , , \ ' scl ' n 



r 2 



x P 



r s d = rj 71 ™- 1 • ■ • fl r sc i = r, 7 ^ n r sol = r„ T i 



(21) 



Computational Efficiency: The time complexity of the greedy 
round robin scheme is heavily based on the computational 
time of the NIntegrate operation in MATHEMATICA. The 
computation time of one NIntegrate operation is denoted 
by r where r = 0.95 sec (on Intel(R) Xeon(R) CPU-X5550 
@2.67 GHz with 24 GB RAM, 64 bit Operating System); this 
is equivalent to the computational complexity of (i) greedy 
scheme and (ii) the first time slot of the greedy round robin 
scheme. Monte-Carlo simulation time required for 100,000 
trials in greedy scheme requires around 150.67 sec which 
is 150 times more than the computational complexity of 
NIntegrate operation. This fact demonstrates the compu- 
tational efficiency of greedy scheme in comparison to Monte- 
Carlo simulations. However, in the second time slot, greedy 
round robin scheme requires K — 1 integrations whereas, for 
the third time slot K — 2 integrations are required. Therefore, 
the computation time at any time slot w < W can be given 
as t + J27=2 T (K — i + 1); w < W, where W denotes 
the total number of time slots. Considering w = W =15 
and the number of rings K = W for greedy round robin, 
the analytical time complexity is around 113 sec which is 
still lower than the Monte-Carlo simulation time required for 
the greedy scheme. Therefore, it can be concluded that the 
numerical implementation of greedy round robin would be 
beneficial for performance assessment in cases with reduced 
number of rings, e.g, when either the number of slots, the 
path-loss exponent, or the cell radius has a small value. 



F. Evaluating the Joint PMF of r sc \ and 9 

Note that, P(r sc \ = r k ) derived for all of the above schedul- 
ing schemes is the marginal PMF of P(r sc \ = r k ,6 = Oi) 
where 9 denotes the angle of the allocated user with respect 
to the serving BS and it is uniformly distributed from to 2ir 
(see Fig. 1). Although the PDF of 9 is continuous we can dis- 
cretize it for analytical consistency and complexity reduction. 
Consider discretizing the range of RV 9 ml uniform angular 
intervals of desired accuracy. Thus P(9 = 9i) = i where 9^ 
denotes any discrete value that the RV 9 can take. Since r sc \ 
and 9 are independent, their joint PMF can be written as: 

P(r sc i = r k ) 



P{r, 



scl 



(22) 



IV. Distribution of Intercell Interference from 
One Cell 

The derivation for the distribution of the ICI from an 
interfering cell I, i.e., fx, (x), depends on the distribution of 
the distance between the allocated user in the interfering cell 
I and the BS of interest, i.e., ff s<!li (r). As mentioned earlier, 
each interfering cell is assumed to have identical conditions in 
a given time slot. Therefore, ff aB i(r) applies to all interfering 
cells and, thus, we will drop the subscript I in the sequel to 
simplify notation. Using the cosine law (see Fig. [TJ, we can 
write: 



,.~2 



D z - 2r sol D cos9 



(23) 



r se i is the distance of the allocated user in the interfering cell 
I from the BS of interest, r sc \ is the distance of the allocated 
user from its serving BS, i.e., (BS /), 9 e [0, 27r] and D = 2R 
since we consider universal frequency reuse with one tier of 
interfering cells. The approach can be extended to any number 
of tiers in a straightforward manner. In order to determine the 
PMF of f sc i where f se i e [D - R, D + R], first of all we 
compute fi t k for given 9 L and r k using d23l as follows: 



,,-2 



^ 2 -2r fc ^cos^ Vr fc ,V0i 



(24) 



where, f^ft denotes the interfering distance from a specified 
polar coordinate (r k ,9i) in the interfering cell to the BS of 
interest located at a distance D (see Fig. H}. In addition, it is 
worth to mention that f,^ are the points at which P(f sc \ = 
fi^ k ) can be defined using d22i i as follows: 



P(r s , 



n, k 



P(r se i = r k ) 



(25) 



The two dimensional data set of f se i, at which P(f se \ = fi >k ) 
is defined, can then be grouped into M segments of any 
arbitrary width A. This can be done by dividing the distance 
between D — R and D + R into M equal segments of width 
and mapping accordingly. Clearly, by adding all the 
probabilities for which f sc i lies in the m th segment we get the 
probability of f sc i = f rn : 



P(n 







E 



P(r so 



Tik) 



(26) 



ri,k£[r„ 



where f m denotes any discrete value that the RV f se i can take. 
Recall X = Kf~^x^ therefore the PDF of X conditioned on 

'Note that A represents the uniform segments of interfering cell whereas 
Afc represents the non-uniform circular regions within a given cell. 
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u< u k f c { lk rlk- l )[F^ lk r p k K^) 



Kr 



-P 



■drfk, (20) 



r sc i can be determined by RV transformation as: 



f. 



A|r„.i 



sel 



(27) 



Averaging over the PMF of f se i, the distribution of the ICI, 
fx{x), from any cell I can be given as: 



fx(x) 



E u 

r m —ri 



P{r S ei = r m ). (28) 



It is important to emphasize that the derivation of the distribu- 
tion of ICI is based on the scheduling decisions of interfering 
cells at a given time slot. Therefore, the parameter r m of the 
ICI distribution varies from one time slot to the other for the 
location based round robin and greedy round robin schemes. 

V. MGF of the Cumulative Intercell Interference 

Computing the distribution of the cumulative ICI Y requires 
the convolution of the PDF of L RVs X u VZ = 1, 2, • • • L, 
which is a tedious task in practice. To avoid convolutions, we 
utilize an MGF based approach to derive the expression for the 
MGF of the cumulative ICI Y. Considering same scheduling 
scheme deployed in each cell, the interference experienced 
from each cell is i.i.d and therefore the MGF of the cumulative 
interference can be calculated as follows: 



M Y (t) = Y[M Xl (t) = (M x (t)) L = {E[e tx ]) J 



(29) 



i=i 



The derivation can be extended in a straightforward manner 
to the case where different cells deploy different scheduling 
scheme^. Looking at the structure of ( l28l . we can write: 



M x (t) = / e tx f x (x)dx 
Jo 



E 



p(r 



sel 



X 







r m —r\ m 
oo 

e tx f x (xriK- l )dx. 



(30) 



The derived expression is generic and applies to any com- 
posite fading distribution. Next, we will present explicit MGF 
expressions for three typically used practical fading models. 
Special Case 1: Rayleigh Fading -(,x ~ Exp(A): In this 
case, the small scale fading coefficient on a given subcarrier 
is considered to be Rayleigh distributed whereas the effect of 
shadowing is not considered. The distribution of interference 
considering a single interfering cell can then be derived as: 



f x (x)= £ K~ 1X 



rie- x ^ R ~ lx P(r Bel = r m ). (31) 



2 The only required change is that in this case, the MGF of ICI received 
from each cell will not be i.i.d, i.e., Mv(t) = nf=i Mx, (t) and the MGF 
of each cell I can be characterized with its corresponding scheduling scheme 
which depends on P(r sc i = rfc). 



Note that (|3H is a Hyper-Exponential distribution with param- 
eter K~ 1 \r!^ l . Thus, using the MGF of the Hyper-Exponential 
distribution, A4y(t) can be derived as follows: 

(. L 
E f> ^~p m - P^ = ?m)\ ■ (32) 

Special Case 2: Generalized-/C Composite Fading x ~ 

/C<3(m s , m c , 0): In wireless channels, shadowing and fading 
across the channel between a user and BS can be jointly mod- 
eled by a composite fading distribution. A closed form com- 
posite fading model, namely Generalized-/C also referred to as 
Gamma-Gamma distribution, has been recently introduced in 
|fT31l which is general enough to model well-known shadowing 
and fading distributions such as log-normal, Nakagami-m, 
Rayleigh etc. Using d28l , fx(x) in this case can be derived 
as follows: 



fx(x) 



, m=fl KfVnT(m c )Y(m s ) 



m c +m s 



(33) 



^m s —trie 



b^xriK- 1 ) P(f se i = f m ), 



where, K„(.) denotes the modified Bessel function of second 
kind with order v, b = 2^/ mc c ™ s ". Performing some algebraic 
manipulations and using 11161 Eq. 6.643/3], the expression for 
Aix(t) can be derived as: 

M x (t) = E p (rsei = ri)e^ ( J 

f m =ri 

& & 2 



AKt 



(34) 

where, W denotes the Whittaker function. Finally, A4y(t) can 
be written as follows: 



M Y (t) 



E ^< 



Tsel 



u ' rn 

AKt 



e sift x 



(35) 



Note that integrating the CDF and MGF of Generalized- 
K, RV which involves Meijer-G and Whittaker functions, 
respectively, in MATHEMATICA and MAPLE can be a bit more 
time consuming. In this context, recently an approximation of 
the Generalized-/C RV has been proposed in ifTBI which is 
discussed below. 

Special Case 3: Gamma Composite Fading x ~ 

Gamma(m s ,m c ): In lfT5l . the authors proposed an accurate 
approximation of the Generalized-/C RV by the more tractable 
Gamma distribution using moment matching method. The 



7 



approximation provides a simplifying model for the composite 
fading in wireless communication systems. Using ( f28l l. fx(x) 
can be written in this case as: 

xr* 3 if- 1 

fx(x) = 22 ,~ ~-B^, m s m. P{r 8e l=f m ). 



is the characteristic function of ( in ring k. The outage 
probability can then be computed by using the classical lemma 
introduced in IfTTl as follows: 



P 



? m =fi #r TO r(m s )m£ 



Im 



(40) 



(36) 



Performing some algebraic manipulations and letting y 



x(^ L — t), Aix{t) can be derived as follows 



where Im((f>z(w)) denotes the imaginary part of 4>z(oj). 
Using d40b . the outage probability can be evaluated using 
any standard mathematical software packages such as 
MATHEMAT I CA and MAPLE. 



M x (t)= £ 



r m —r 



x r(m s ) (fiK- 1 - m c t) ms -/o 



e a y" is X (iy Evaluation of Ergodic Network Capacity: Another 
important performance evaluation parameter is the network 



E ^=r m ) - 



(37) 



Finally, .M^i) can be written as follows: 

f ^ ( K- 1 ^ 

M r{t) = E = f - * 

\f m =f 1 



K 1 fm — m c t 



(38) 



VI. Evaluation of Important Network 
Performance Metrics 

In this section, we demonstrate the significance of the 
derived MGF expressions of the cumulative ICI in quantifying 
important network performance metrics such as the outage 
probability P ou t, ergodic capacity C and average fairness T 
among users numerically. 

Evaluation of Outage Probability: The outage probability 
is typically defined as the probability of the instantaneous 
interference-to-signal-ratio to exceed a certain threshold. In 
order to evaluate P out numerically, we use the MGF based 
technique introduced in |[T7l for interference limited systems. 
Firstly, we define a new RV, Z — qYlt=i X — Xo = qY — Xq, 
where q is the outage threshold and Xq is the corresponding 
signal power of the scheduled user in the central cell. 
An outage event occurs when p{Z > 0), i.e., when the 
interference exceeds the corresponding signal power. This 
decision problem is solved in IfTTl by combining the 
characteristic function of Z and residue theorem. The 
characteristic function of Z is defined as cf)z(ju) = E[e jZ "]. 
Considering interference Y and signal power Xq to be 
independent, the expression for <f>z(juj) can be given as, 
<f>z{ju) = <fc{jqu)<j)x Q {-ju); where 0y(qjw) can be given 
by (|32| |. ( l35l l. and ( f38l ) for different fading models. In general, 
the characteristic function of Xq can be calculated as |fT8l : 

/>oo 

4> Xo (w) = E(e^°) = / ei" x °f Xo (x )dx 

}°oo (39) 

= ju e^°F Xo (xo)dxo, 
Jo 

where the detailed expressions of Fx (xo) for all opportunistic 
scheduling schemes are provided in iflOl . Compact closed 
form expressions of <fix (jw) are available in the literature for 
e.g. |fl9l Eq. 19]. For non-opportunistic scheduling schemes 

^Xo(jw) = Erh=n < / , Ck fc (jw)i 3 (r sc i = r k ), where <f>c\ rk (ju) 



lo£ 



1 



ergodic capacity C, i.e., C = ^ 

L V 2^i=i A 
Usually, the computation of C requires (L + l)-fold numerical 

integrations. To avoid this, we utilize the efficient lemma 

derived in ll20l with a slight modification to take thermal 

noise into account as follows: 



E 



In 1 



Xo 



M Y (t) -M Xo , Y {t) 



(41) 



t 



*dt, 



where, M Y (t) = E[e-*^i=i Xi ] and Mx Q .y{t) = 
E [ e -t(x +Ef=iX,)] = E[e-*( Xo+Y )]. Note that this is the 
definition of MGF as defined in l20l which is not the same as 
our definition. Thus, we can use My it) from (|32| |. ( f35l > and 
( l38t directly with a sign change of jw. In addition, fiTI ) can 
be solved efficiently by expressing it in terms of the weights 
and abscissas of a Laguerre orthogonal polynomial l20l as 



E 



In 1 



+ 



Re, where ^ e and a € are the sample points and the weight 
factors of the Laguerre polynomial, tabulated in II 1 611 . and 
Re is the remainder. The MGF of Xo can be calculated as 
explained in (|39l . 

Evaluation of Average Fairness: In order to quantify the 
degree of fairness among different scheduling schemes, we 
use the notion developed in lf2Tl . The average fairness of 
a scheduling scheme with U users can be given as, T = 



— J2i=i Pi log 10 77 » wnere Pi is the proportion of resources 
allocated to a user i or the access probability of user i, 
A system is strictly fair if each user has equal probability 
to access the channel and in such case the average fairness 
becomes one. The other extreme occurs when the channel 
access is dominated by a single user; in such case, the average 
fairness reduces to zero. The average fairness can be easily 
computed using our derived results as 



T 



K 

E 

k=l 



Pfr« 



rk) 



lo gio p ( r sci = rk) - log 10 u fe 



log 10 E7 



(42) 



where u k denotes the number of users in a ring k. 

VII. Numerical and Simulation Results 

In this section, we first define the system parameters and 
discuss the Monte-Carlo simulation setup which is required to 
demonstrate the accuracy of the derived expressions. We then 
address some important insights and study the performance 
trends of the different scheduling schemes. 
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A. Parameter Settings and Simulation Setup: 

The radius R of the cell is set to 500m and the cell is 
decomposed into non-uniform circular regions of width Afc. 
The path loss variation within each circular region is set to 
K = 2dB. 

For each Monte-Carlo trial, we generate U uniformly 
distributed users in a circular cell of radius R. Each user 
has instantaneous SNR given by (Q]i and short term average 
selected SNR based on its ring location (7^). We allocate a 
user with maximum instantaneous SNR in the greedy scheme, 
a user with maximum normalized SNR in the proportional 
fair scheme and a user arbitrarily for the round robin scheme. 
For location based round robin we select a user randomly 
from the w th ring in a time slot w whereas we select a 
user with maximum SNR without considering the users of the 
previously allocated w — 1 rings for the greedy round robin 
scheme. Next, we calculate the distances of the selected users, 
i.e., r se i, from their serving BS and compute the distances to 
the BS of interest, i.e., r sc \ for all scheduling schemes. The 
process repeats for large number of Monte-Carlo simulations. 
The distance data is then analyzed by creating a histogram 
of non-uniform and uniform bin width in Fig. [2] and Fig. [3] 
respectively. 
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Fig. 2. PMF of the distance of the allocated users in a given cell (i.e., PMF 
of r se i) for proportional fair, greedy, and round robin scheduling schemes 
with path loss exponent /3 = 2.6, U = 50, C=60 dB, P ma x=lW, (t 2 =-174 
dBm/Hz, and Number of Monte-Carlo simulations =100,000. 



B. Results and Discussions: 

Fig. |2] depicts the PMF of the location of the scheduled 
user in a given cell based on the proportional fair, greedy 
and round robin scheduling schemes. Since the proportional 
fair scheme exhibits some fairness among users in a cell, the 
PMF of the allocated user locations is expected to be more 
flat compared to the greedy scheme. Since the cell edge has 
more users due to the large area and each user has equal 
probability to be allocated on a given subcarrier, therefore 
the round robin scheme exhibits high probability at the cell- 
edge. In order to get an integer number of users within a 
ring, we perform rounding in the analysis, i.e., we consider 
zero active users in the rings where Uk < 0.5. In Monte- 
Carlo simulations, we consider the probability of allocating 
a user in these rings to be zero which can also be verified 
from Fig. |2] The impact of rounding can be observed more 
in the cell center than at the cell-edge, due to narrow circular 
regions and in turn low number of users in each ring in the 
cell center. This rounding effect is therefore more visible in 
the PMF of the greedy scheme which exist mainly in the cell 
center as observed in Fig. 2(b). However, as this mismatch 
lies near the cell center, it does not have much impact on the 
derived performance results which include distribution of the 
ICI, ergodic capacity and outage probability. 

It is important to note that the numerical results for the 
derived PMF in Fig. [2] nearly coincide with the exhaustive 
Monte-Carlo simulation results with a small number of rings 
K = 10 and k = 2dB. Moreover, it can also be noticed 
that the width of the circular regions tend to increase from 
cell center to cell edge which is due to the exponentially 
decaying path loss as mentioned in Section II. The number 
of required rings is expected to decrease by reducing f3 and 
increasing the amount of power decay within each circular 
region and vice versa. Another important point to explain 
with reference to Fig. [2] is that with the increase in the 
number of competing users on a given subcarrier, the PMF of 



opportunistic scheduling schemes will get more skewed toward 
the cell center which is due to the fact that the higher the 
number of users in the cell center, the higher is the probability 
of allocating a subcarrier in the cell center. 

In Fig. [3] the PMF of the distance between the allocated user 
in interfering cell I and the BS of interest, i.e., P(r sc \ = f m ), 
is presented. Numerical results are found to be in close agree- 
ment with the Monte-Carlo simulation results and M = 20, 
i.e., A=50 m. For the opportunistic scheduling schemes, it is 
likely that a user close to its serving BS can get a subcarrier, 
thus, the PMF of the distance of allocated interfering users 
is expected to have high density in the middle. However, the 
slight descend in the central region in Fig.[3]is due to ignoring 
users that lie within the rings where the average number of 
users is less than half. Moreover, we can observe that the round 
robin scheduler is highly vulnerable to interference compared 
to the other schemes as high interference is expected to come 
from the cell edge users in the interfering cells. On the other 
hand, the greedy scheduler is expected to have allocations 
near the cell center and, thus, leads to less interference from 
neighboring cells. The proportional fair scheme lies in between 
the two extremes. Note that as M increases the analytical 
accuracy is expected to increase even further with a trade-off 
in terms of computational complexity. An acceptable range of 
A with negligible complexity varies from 1 to 100 m. 

Fig. [4] illustrates the CDF of the ICI considering different 
number of interfering cells and path loss exponents (3 for the 
greedy scheduling scheme. With the increase in the number 
of interferers, the interference level increases. Moreover, as 
j3 increases, the signal degrades rapidly and thus interference 
level is reduced considerably. At this point, it is important to 
mention that in this paper we derive and utilize the MGF of 
the cumulative ICI rather than the CDF of the cumulative ICI 
in order to evaluate important network performance metrics. 
Therefore, the analytical part of the provided figure of the CDF 
of the cumulative ICI is plotted using a technique mentioned 
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Fig. 3. PMF of the distance at which the users in the interfering cells are 
allocated (i.e., PMF of r se i) for proportional fair, greedy and round robin 
scheduling schemes with path loss exponent /3 = 2.6, U = 50, X = 180, 
X ~ Gamma(3/2,2/3), C=60 dB, P max =lW, <r 2 =-174 dBm/Hz, A=50 
m, and Number of Monte-Carlo simulations =100,000. 



Fig. 4. Impact of different number of interferers L and various path loss 
exponents on the CDF of cumulative ICI considering greedy scheduling 
scheme, U = 50, X = 180, \ ~ Gamma(3/2, 2/3), C=60 dB, P max =lW, 
cr 2 =-174 dBm/Hz, and Number of Monte-Carlo simulations =100,000. 



in 1221 to convert MGF into CDF numerically. 

Fig. [5] investigates the effect of increasing the number 
of competing users on a given sub-carrier considering all 
scheduling schemes. It can be observed that the increase in the 
number of users enhances the performance of the opportunistic 
scheduling schemes due to additional multiuser diversity gains. 
The greedy scheme achieves the best performance whereas 
the round robin scheme achieves the worst performance. As 
expected, the proportional fair scheme lies in between the two 
extremes. The average capacity of location based round robin 
over W = K time slots has been shown to be better than 
the conventional round robin scheme. The average capacity 
results of the greedy round robin scheme is presented for 
W = 3 and W = 6. Clearly, for W = 1, the scheme is 
equivalent to the greedy scheme; however, with the increase 
of time slots, performance degradation takes place due to the 
reduction of multiuser diversity caused by ignoring the users 
from previously allocated rings. 

Fig. [6] quantifies the average resource fairness of all pre- 
sented scheduling schemes. As expected, round robin is a 
strictly fair scheme. The proportional fair scheme possesses the 
ability to enhance the network throughput compared to round 
robin scheduling while providing a high degree of fairness. 
The greedy scheme is observed to be the most unfair scheme. 
Considering K time slots, the average fairness of the location 
based round robin scheme is investigated and found to be very 
close to the round robin scheme, however, with degradation 
in performance as can be observed in Fig. [5] For the greedy 
round robin scheme, we plotted the fairness metric considering 
W = 3 and W = 6; it is shown that as the number of time 
slots increases, the fairness improves with a trade-off price in 
terms of ergodic capacity. 

In Fig. [7] we evaluate the network outage probability as 
a function of the outage threshold; q = (Z + Xq)/Y for 
(i) U — 50 users; (ii) U = 100 users. The numerical and 
simulation results are nearly identical for most cases. The 
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Fig. 5. Network ergodic capacity for different number of users considering 
different scheduling schemes with ft = 2.6, x ~ Gamma(3/2, 2/3), C=60 
dB, P max =lW, and o- 2 =-174 dBm/Hz. 



higher the outage threshold for a given signal and interfer- 
ence power, the greater outage is expected. Moreover, for 
larger number of users the outage probability is observed to 
reduce for all opportunistic scheduling schemes except the 
round robin scheme. Since increasing the number of users 
on a given subcarrier in non-opportunistic schemes does not 
directly affect the access probability of a ring k, therefore 
its impact on the ICI is almost negligible. This fact can also 
be verified from Fig. [5] Finally, in Fig. [8] we evaluate the 
network ergodic capacity as a function of the fading severity 
parameter and average power of Gamma fading interference 
channels for different scheduling schemes. Firstly, it can be 
observed that increasing the average power of the interference 
channel which is given by ft = m c m s for a given fading 
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Fig. 6. Average system fairness for different number of users considering 
different scheduling schemes with (3 = 2.6, 0=60 dB, P m ax= 1 W, and <r 2 =- 
174 dBm/Hz. 



Fig. 8. Ergodic network capacity of different scheduling schemes as 
a function of different parameters of the interference statistics \ ~ 
Gamma(m,,m c ) with U = 50 users, ,8=2.6, and ct 2 =-174 dBm/Hz. 
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Fig. 7. Impact of different scheduling schemes on the network outage 
probability for different number of users and various outage thresholds, (3=2.6, 
X ~ Gamma(3/2, 2/3), (7=60 dB, P max =lW, and <r 2 =-174 dBm/Hz. 

severity parameter m s , the capacity degrades significantly for 
all schemes. Moreover, it is also shown that increasing the 
fading severity m s while keeping the average power = 3 
fixed has minimal impact on the system capacity. Therefore, 
the lower average power of interference channel S7, the better 
is the overall system performance. 

The small gap between the analytical and simulation results 
is mainly due to assuming in the analytical derivations that 
users located within a ring are at the boundary of the ring. 
This gap can be further reduced in Fig. [2] by increasing the 
accuracy of the approximation, i.e., by reducing the path loss 
decay n within each circular region which in turn increases 
the number of rings. An acceptable range of k with reasonable 
complexity varies from 0.5 to 3 dB such that r < R. 



VIII. Conclusion 

We proposed a novel approach to model the uplink ICI 
considering various scheduling schemes and composite fading 
channel models. The proposed approach is not dependent on 
a particular shadowing and fading statistics, hence, extensions 
to different models is possible. The provided numerical results 
help in gaining insights into the behavior of ICI considering 
different scheduling schemes and composite fading models. 
Moreover, they provide quantitative assessment of the relative 
performance of various scheduling schemes which is important 
for network design and assessment. The proposed approach 
can be extended to typical power control schemes and down- 
link scenarios as discussed in [10]. 
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